Entropy fluctuations for directed polymers in 2+1 dimensions 



O 

o 
o 

(N 
-i— > 
O 

O 



I 

C/3 



I 

C 

o 
o 



On 
(N 
O 

o 
o 

c3 



c 

o 
o 



X 



Xiao-Hong Wang 1 ' 2 , Shlomo Havlin 1 and Moshe Schwartz 3 
1 Minerva Center and Department of Physics, Bar-Ban University, Ramat Gan 52900, Israel 
2 Department of Thermal Science and Energy Engineering, University of Science and Technology of China, Hefei, Anhui 

230026, P. R. China 

' School of Physics and Astronomy, Raymond and Beversly Sackler Faculty of Exact Science, Tel-Aviv University, Tel-Aviv 

University, Tel-Aviv 69978, Israel 
(February 1, 2008) 



We find numerically that the sample to sample fluctuation 
of the entropy, AS, is a tool more sensitive in distinguishing 
how from high temperature behaviors, than the correspond- 
ing fluctuation in the free energy. In 1+1 dimensions we find 
a single phase for all temperatures since (AS) 2 is always ex- 
tensive. In 2+1 dimensions we find a behavior may look at 
first sight as a transition from a low temperature phase where 
(AS) 2 is extensive to a high temperature phase where it is 
subextensive. This is observed in spite of the relatively large 
system we use. The observed behavior is explained not as a 
phase transition but as a strong crossover behavior. We use 
an analytical agreement to obtain (AS) 2 for high temperature 
and find that while it is always extensive it is also extremly 
small and the leading extensive part decays very fast to zero 
as temperature is increased. 
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The problem of directed polymers has attracted much 
interest in recent years. It is relevant to many fields 
ranging from surface growth phenomena and spin glasses 
to flux lines in high-T c superconductors [1-4]. It is well 
known that the problem of directed polymers in a random 
medium is equivalent to the KPZ equation that describes 
surface growth [5-7]. Much is known about the KPZ sys- 
tem, in particular in 1+1 dimensions. The KPZ equation 
provides the exact dynamical exponent for directed poly- 
mers in 1+1 dimensions [5-9]. The situation for higher 
dimensions is more complex. Traditional approximation 
schemes like dynamical renormalization-group methods 
fail to produce the exponents obtained by simulations 
[1,2, 10-15]. A self-consistent expansion of the correla- 
tion function introduced later [16] yields results compat- 
ible with simulations for 2+1 dimensions. Above 2+1 di- 
mensions, the behavior of directed polymers has a phase 
transition when temperature is raised (in the directed 
polymer problem or the level of noise in the KPZ sys- 
tem) and the system goes over from a strong coupling 
to a perturbative weak coupling behavior [17-20]. All 
the field theoretical treatments agree that at 2+1 dimen- 
sions, namely at the lower critical dimension itself, no 
transition should occur and the strong coupling behav- 
ior exists at zero and any finite temperature [16,19-21]. 
Some authors have claimed in the past on the basis of nu- 



merical simulations, to have obtained a phase transition 
in 2+1 dimensional systems [22-28]. The system studied 
were, however, relatively small and all those claims were 
not pursued eventually. In this paper, we present numer- 
ical results obtained for much longer systems in 1+1 and 
2+1 dimensions. We study the free energy fluctuations 
that is usually studied in literature but in addition we ob- 
tain numerically also the sample to sample fluctuations 
of the entropy. This is a quantity, that is less common in 
the literature and was introduced first, to the best of our 
knowledge, by Fisher and Huse [29]. The fluctuations in 
the entropy are more pronounced than those of the free 
energy and in fact, we expect it to be extensive in the size 
of the polymer in 1+1 and 2+1 dimensions. In higher di- 
mensions we expect a transition from a low temperature 
extensive phase to a high temperature subextensive one. 
Namely defining the quantity X(t) — (AS) 2 /t, its infi- 
nite volume limits, behaves like an order parameter. It is 
zero at high temperatures and is of order 1 at low tem- 
peratures. The numerical results concerning X(t) in 1+1 
dimensions are those being expected. We find that A is of 
order 1 over all the temperature region. In 2+1 dimen- 
sions, the numerical results seem to indicate an unex- 
pected transition from a low temperature phase in which 
A is of order 1 into a high temperature regime where it is 
approaching zero. Although the transition temperature 
seems to be defined very sharply and although X(t) can 
be fitted at high temperatures by X(t) cx Int/t, we claim 
that it is just a crossover phenomenon and the reasons 
for it occuring in such a spectacular way will be discussed 
later. 

Consider a directed polymer on a hypcrpyramid lat- 
tice structure with the random energy assigned on each 
bond. The partition function G(R, t) for directed poly- 
mers starting from (0, 0) and ending at (R, t) is defined by 

G(R, t) = J2 e.~ Ec / T where Ec is the sum of the energy 

c 

on the path C and T is the temperature. For simplic- 
ity, we demonstrate our calculations using the transfer- 
martrix method for the case of 1+1 dimensions. A sim- 
ilar formalism has been used for 2+1 dimensions. The 
iteration relation for the partition function G(R, t) is 



G{R, t + l) = G(R - 1, i)e~ £i/T + G(R + 1, t)e~ er 
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in which, ei and e r are the energy assigned to the 
left and right bonds of the point (R,t). The free en- 
ergy F(t) is given by F(t) = -TlnG(t), where G(t) = 

£ G(R, t) is the total partition function. The free en- 
R 

ergy fluctuation AF = (F 2 — F ) 1 / 2 has been com- 
monly studied (A is the ensemble average of the quan- 
tity A). We can also define the internal energy, (E) = 

£ £ F c e- EcjT I £ £ e- Ec / T . The internal energy flue 
ft C ft c 

tuation (AE) T = ((E) 2 - JE) 2 ) 1 / 2 is also an inter- 
esting quantity. In order to obtain the iteration rela- 
tion for the internal energy (E), we define E(R,t) = 
£ E C (R,t)e' Ec ^ t)/T /G(t). It is clear that (E) = 

C(R,t) 

Y^E(R,t). The iteration relation for E(R,t) is: 
R 

E(R, t + l) = [er^l T E(R - 1, t)G(t) + e -^/ T E(R + 1, t) 
G(t) + e ie - e '/ T G(R - 1, t) + e r e- e "/ T G(R + 1, t)]/G(t + 1) 

(2) 

To understand the difference between the free energy 
and internal energy fluctuations, we use the entropy 
S = ((E) — F)/T. We present systematic numerical sim- 
ulations for the entropy fluctuation AS = (S 2 — S ) 1//2 
based upon equations (1) and (2) with the initial con- 
ditions G(R,0) = 5 R , and E(R,0) = 0. The random 
energy assigned on the bond is assumed to be uniformly 
distributed in the interval (-0.5,0.5) and uncorrelated in 
space and time. We use length up to t = 2000 (d = 1 + 1) 
and t = 1000 (d = 2 + 1) [30]. Six thousands configu- 
rations for 1+1 dimensions and four thousands config- 
urations for 2+1 dimensions were collected to take the 
ensemble average. 

In 1+1 dimensions, the numerical results clearly show 
that the entropy fluctuation has the behavior (AS) 2 oc t 
for any temperature (see Fig. 1) [30]. It is expected 
that the entropy fluctuation will tend to zero at the lim- 
its of zero and infinite temperature. Indeed, the slope of 
(AS) 2 /t (Fig. 3(a)) is about 1 for T < T P , i.e. \(t) ~ T 
and about -4 for T > T P , i.e. \(t) ~ T~ 4 . As a 
result, the free energy fluctuation and internal energy 
fluctuation will be the same at the two limits. About 
T = Tp = 0.2, the entropy fluctuation reaches a maxi- 
mum. There is no evidence for a phase transition in 1+1 
dimensions. 

The picture of directed polymers in 2+1 dimensions 
is more complicated than that of 1+1 dimensions (see 
Fig. 2). Similar to 1+1 dimensions, the entropy fluc- 
tuation tends to zero at the two limits of zero and infi- 
nite temperature and there is a peak at Tp = 0.11. For 
low temperatures T < Tp = 0.11, we see that the en- 
tropy fluctuation (AS*) 2 tends to t and (AS) 2 ~ T as 
in 1+1 dimensions (see Fig. 2(a) and Fig. 3(b)). How- 
ever, for T > Tp, we find that the increase of entropy 



fluctuation as a function of t becomes slower and slower 
as temperature is increased (see Fig. 2(b)). At the very 
high temperatures, e.g. T — 10.0, the entropy fluctu- 
ation is proportional to Int for large t. It seems that 
(AS) 2 /t will tend to the nonzero values only in the re- 
gion < T < Tp. The fluctuations of the free energy AF 
are correlated with the fluctuations of the entropy AS. 
The high temperature behavior of AF gives for large t a 
logarithmic dependence on t while for lower temperatures 
we see a crossover from logarithmic behaviour at small t 
to AF ~ t - 2 for larger i's. We have not seen however, a 
sharp obvious temperature where this happens. In Fig. 
3b we see a sharp transition from a low temperature re- 
gion T < Tp, where X(t) is almost t independent to a 
high temperature region where \(t) is a decreasing func- 
tion of t. The transition temperature is very close to the 
point where A(T, t) is maximal as a function of T for all 
t. 

The observed behavior may be explained as follows. 
The infinite system is characterized by a correlation 
length £. As long as the size of the system is adequately 
described in terms of the linear theory of deposition pro- 
posed by Edwards and Wilkinson [31]. It is easy to 
show that within the Edwards- Wilkinson theory (AS) 2 is 
proposional, for a finite system to Int, as obtained by us 
for high temperatures. For t longer than £ the nonlincar- 
ities become important and (AS) 2 should be extensive in 
t. How is this related to the temperature dependence we 
find? The correlation length £ has a strong dependence 
on temperature, £(T) is proposional to exp[(T /T o ) ] with 
6 = 3 according to Fisher and Husc [29] and 6 = 2 accord- 
ing to Kim, Bray and Moore [18]. In any case the tem- 
perature dependence of £t is so strong that a relatively 
small increase in temperature may result in an increase 
of £t by orders of magnitude increasing it from values 
below the minimal t we are using (t = 50) to well above 
the maximal value (t = 1000). Indeed, a more careful 
examination of the data is consistent with the above ex- 
planation. In the temperature region Tp < T < 2Tp, 
all the lines but the one corrcspondings to the smallest 
t merge. This suggest that the asymptotic value of X(t) 
has already been reached, so that the asymptotic value 
is of order 1. Thus above Tp the asymptotic value is 
still finite. This may explain the former numerical re- 
sults that claimed a phase transition in 2+1 dimensions 
[22-26]. The strong dependence of the correlation length 
on temperature suggest that increasing the size of the 
system does not really undergo a phase transition. 

It is not difficult to obtain the asymptotic (infinite t) 
high temperature form of X(T). As expected from the 
Edwards- Wilkinson model X(t) decrease as Int/t. This is 
true as long as t > £t but as t becomes of the order of £tA 
must tend to a constant independent on t. Therefore, 

A(T) cc ^ = (3) 
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Thus we see that although A is not zero it decreases ex- 
tremely fast with temperature. 
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FIG. 3. Plot of entropy fluctuation per unit length of the polymer as a function of temperature (a) for the different time 
t = 50,500, 1000,2000 in 1 + 1 dimensions, (b) for the different time t = 50, 100,300,600, 1000 (from top to bottom) in 2+1 
dimensions 
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